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Abstract 

We examine some nontrivial consequences that emerge from interpreting a 
■ position-dependent mass (PDM) driven Duffing oscillator in the presence of a 

quartic potential. The propagation dynamics is studied numerically and sensi- 
' tivity to the PDM-index is noted. Remarkable transitions from a limit cycle to 

& • chaos through period doubling and from a chaotic to a regular motion through 

intermediate periodic and chaotic routes are demonstrated. 
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I/* - ) ■ While position-dependent (effective) mass (PDM) quantum mechanical systems 

have repeatedly received attention in many areas of physics (such as for instance, 
in the problems of compositionally graded crystals [T], quantum dots [2], nuclei [3], 
quantum liquids jl], metal clusters [5] etc; see also [HI El [91 [101 EEH H21 EE] for 
theoretical developments and references therein), interest in classical problems having 
a PDM is relatively recent and rapidly developing subject [M [T51 IM [TT1 [TBI [P31 12U1 
EU [22J [23]. The simplest case of a PDM classical oscillator has been approached by 
several authors [El EH HH1 US EQ] • In what follows we shall demonstrate that a PDM 
ascribed Duffing oscillator provides an attractive possibility of defining a dynamical 
system that exhibits various features of bifurcations, chaos and regular motions. To 
be specific we will focus on the following equation of motion which emerges from a 
Lagrangian of the form [21] 

L = ^ ( TT^ )(i;2 " a;oV) ' eeM (1) 
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(where an overhead dot indicates a time derivative) namely, 

(1 + ^x 2 )x - £xx 2 + u$x = 0. (2) 

Equation (2) represents a special type of nonlinear oscillator which can 
be reduced to a first order form by effecting a change in the variable 
x | I24] . The latter equation can be immediately solved yielding periodic 
solutions. Hence the absolute regularity of the motion does not depend 
on the specific value of £. 

Classically when the mass is position dependent, Newton's equation of motion 
gets modified to 

m(x)x + m\x)x 2 = (3) 

(where the prime indicates a spatial derivative) in the absence of any external force 
term. When compared with fl2]) the following profile of the mass function comes out 

mix) := — . (4) 

on ignoring the presence of the harmonic term u 2 x to effect such a comparison. In 
d3J we have scaled the constant mass to unity. Notice that it is also possible to deal 
with other mass functions [TJ)J [2D] but for concreteness we focus on (j3J) in the rest of 
this paper. 

In [20] an attempt was made to construct the underlying Lagrangian wherein the 
PDM system (jHJ) is acted upon by a force F that could be dependent on position x, 
velocity x and time t i.e. 

dp 

F(x, x, t) := — = m\x)x 2 + m(x)x (5) 

where p = m(x)x is the linear momentum. They have found that, in addition to the 
kinetic energy T, a reacting thrust R is at play as well 

T:=-m(x)x 2 , R(x, x, t) := — -m'{x)x 2 . (6) 



where R is essentially a non-inertial force. Assuming that F can be split as F - 
F{x) + R(x,x,t), where F(x) is controlled by a scalar potential function F = — |g 
the Lagrangian form of (J5J obeys 



Corresponding to the above L which is L = |m(x)x 2 — V(x), the Hamiltonian reads 
H = 2m{x) ^ V( x )- These are of standard text-book formal with m = m(x). 



1 Quantum mechanically such forms are evidently non-Hermitian 
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We stress that the non-potential force term R has to depend on the velocities 
for otherwise they can be identified with the potential force F(x). When written 
explicitly ([7]) looks 

d dT _ dT _ _dV ^ 
dt dx dx dx 
As an aside, the time rate of change of T can be easily worked out as 

dT d .dT .. ,dV dT 

-r ^ -r(^rx) + (- R)x + — 9 

dt dVdx 1 y dx ' dt w 

from which, for a scleronomic system it at once follows that the time rate of energy 
E{= T + V) is given by 

dE ~ 1 , , 1 fxx 3 , , 

= Rx = — m'x 3 = . 10 

dt 2 2(1+^2)^/1 + ^2 

fTTOT) speaks of the power of the non-potential force. This result is new for PDM 
systems and of interest. If the power is negative we encounter dissipative systems. 

The purpose of this communication is to develop a mathematical framework that 
enables us to study the evolution of the PDM system flS]) in the presence of an 
external periodic (non-autonomous) force with an additional damping term moving 
in a quartic potential. Towards this end we focus on an extended PDM equation of 
motion 

m{x)'x + m\x)x 2 + oJ^x + Xx 3 + ax = f cosut (11) 

which reduces to ([3]) in the absence of all the force terms. 

Some remarks are in order [261 EH EH] • For the constant mass case i.e. m(x) = 1, 
ffTTT) goes over to a forced, damped Duffing oscillator which because of the presence 
of a double-well potential mimicks a magneto-elastic mechanical system. The latter 
is concerned with a beam placed vertically between two magnets with a fixed top 
end and free to swing at the bottom end. As soon as a velocity is enforced the beam 
begins to oscillate eventually coming to rest at an equilibrium point. However the 
situation changes when a periodic force is applied : stable fixed points or stable fixed 
angles no longer occur. In (iTTj) A is the governing parameter along with the periodic 
force f cosut in the presence of a viscous drag of coupling strength a. 

Taking m(x) as in (j3J) we obtain from fill I) the coupled set of equations 

x = y 

l + £x 

z = w 



y = 1 ; Cm2 + a/1 + ix 2 [jcosz - u x - Xx 3 - ay] (12) 



The remarkable behavior of the dynamical system (fT2|) can be understood by 
examining the interplay between the amplitude / of the periodic forcing term and 
the PDM parameter £. Various studies of the corresponding constant-mass case, i.e. 
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when £ = 0, have been made in the literature (see specifically [29]) and the nonlinear 
behavior demonstrated including the oscillation modes and the nonlinear resonances, 
both theoretically and experimentally. For the PDM case, we fix the parameter values 
to be oj = 1.0, Wq = 0.25, a = 0.2 and A = 1.0 as is standard [23]. Further, we will 
always assume £ > 0. We plot in Figure 1 the phase diagrams for different values 
of £ including the constant-mass case of £ = 0. While in the latter situation we 
encounter limit-cycle oscillation, there is a drastic change in the dynamical behaviors 
as £ is increased. For instance, period-four oscillations are observed for £ = 0.2 which 
eventually give way to a chaotic behaviour both for £ = 0.4 and £ = 0.6 values. 

Bifurcation diagram of the system (1121) for x with respect to / by taking £ = 0.5 
are presented in Figure 2. As soon as periodic force is applied, limit cycle oscillations 
gather in the range < / < 2. But with the increase of the amplitude of the forcing 
term, period-two oscillations set in and survive briefly up to / = 4. Further stepping 
up of / produces period-four and period-eight oscillations ultimately leading to a 
chaotic behaviour. But such a regime is short-lived because with the /-value going 
beyond 7.5, a period-halving bifurcation appears yielding period-three oscillations for 
/>8. 

In Figure 3 we demonstrate the bifurcation sequence for x but with respect to 
£ by taking / = 5.0. As is evident, the limit cycle turns to period-four oscillations 
progressing to period-eight and so on running into a chaotic behaviour that lasts for 
a while before a regular motion takes over. The latter then again yields to a chaotic 
dynamics and the entire motion subsequently becomes a regular one around £ = 1.9. 
On the other hand, by carrying out a bifurcation analysis for x with respect to £ 
by taking / = 8.0 we find from Figure 4 that the constant-mass oscillator reveals a 
chaotic state. But we also notice that a subtle interplay between the parameters / 
and £ produces complicated dynamics from a chaotic phase to periodic oscillations 
back again to a chaotic character and finally settling into a regular bahaviour. 

To summarize, we have demonstrated in this communication that given any value 
of /, if it is confronted with the PDM parameter £, remarkable phase transitions occur 
such as, for instance, from a limit cycle mode to a chaotic regime through a period 
doubling intermediate phase pointing to the sensitivity of the system ( II IB due 
to the interplay between / and £. The essential distinction between the 
constant mass and the variable mass case rests in the fact that the presence 
of the parameter £ not only enhances the rapidity of such transitions but 
also initiates complicated nature of bifurcations, of course for a non-zero 
value of /. Conversely, we also encounter transitions such as from a chaotic phase 
— > periodic orbit — > period doubling — > chaos — > regular motion. 

We thank the anonymous reviewers for their constructive suggestions. 
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Figure 1: Phase diagram of equation (TT2"j) with / = 5.0, u = 1.0, 0Jq = 0.25, a = 
0.2, A = 1.0 for (a) constant mass (f = 0)(b) PDM with £ = 0.2 (c) PDM with 
£ = 0.4 (d) PDM with f = 0.6 




Figure 2: Bifurcation diagram of equation ffl2|) with respect to / for £ = 0.5, u; = 
1.0, = 0.25, a = 0.2, A = 1.0 
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Figure 3: Bifurcation diagram of equation (|12p with respect to £ for / = 5.0, u; 



1.0, = 0.25, a = 0.2, A = 1.0 
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Figure 4: Bifurcation diagram of equation (fl2~]) with respect to £ for / = 8.0, u 



1.0, ul = 0.25, a = 0.2, A = 1.0 
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